Hot QCD equations of state and response functions for quark-gluon plasma 
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We study the response functions (chromo-electric susceptibilities) of quark-gluon plasma as a 
function of temperature in the presence of interactions. We consider two equations of state for 
hot QCD. The first one is fully perturbative, of 0{g'') EOS and, and the second one which is 
0[g^ \n{l/ g) + S\, incorporates some non-perturbative effects. Following a recent work (Physical 
Review C 76, 054909(2007)), the interaction effects contained in the EOS are encapsulated in 
terms of effective chemical potentials(/i) in the equilibrium distribution functions for the partons.By 
using them in another recent formulation of the response functions! arXiv : 0707 . 3697 ) , we determine 
explicitly the chromo-electric susceptibilities for QCD plasma. We find that it shows large deviations 
from the ideal behavior. We further study the modification in the heavy quark potential due to 
the medium effects. In particular, we determine the temperature dependence of the screening 
lengths by fixing the effective coupling constant Q which appears in the transport equation by 
comparing the screening in the present formalism with exact lattice QCD results. Finally, we 
study the dissociation phenomena of heavy quarkonium states such as cc and bb, and determine the 
dissociation temperatures. Our results are in good agreement with recent lattice results. 

Keywords: Response function; non-Abelian permittivity; Quark-Gluon plasma; hot QCD 
equation of state; equilibrium distribution function; chemical potential; RHIC. 
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I. INTRODUCTION 

It is expected that at high temperatures (T ^ 150 — 
2QQMeV) and high densities (p ^ IQGev / fw?) nu- 
clear matter undergoes a deconfinement transition to the 
quark-gluonic phase. This phase is under intense inves- 
tigation in heavy ion collisions, and already, interest- 
ing results have been reported by Relativistic Heavy Ion 
Collider(RHIC) experiments [l|. As an important de- 
velopment, flow measurements suggest that close to 
the transition temperature Tc, the quark-gluon plasma 
(QGP) phase is strongly interacting — showing an al- 
most perfect liquid behavior, with very low viscosity to 
entropy ratio — rather than showing a behavior close to 
that of an ideal gas. See Ref. Q for a comprehensive re- 
view of experimental observations from RHIC, and Ref. 
[IIjSHjSH for other recent experimental results. On the 
other hand; lattice computations [1, Q also suggest that 
QGP is strongly interacting even at T = 2Tc. This find- 
ing has been reproduced by a number of other theoret- 
ical studies — by employing AdS/CFT correspondence 
in the strongly interacting regime of QCD [3, by molec- 
ular dynamical simulations for classical strongly coupled 
systems and by model calculations with Au-Au data 
from RHIC [H 

If this be the case, as it indeed appears to be, then 
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the plasma interactions would be largely in the non- 
perturbative regime; in this regime, few analytic tech- 
niques are available for a robust theoretical analysis. Ef- 
fective interaction approaches are needed. In this direc- 
tion, considerable work has already been done and we 
refer the reader to Ref. [H, [13, [H, 111, S [HI, [H, ^ for 
some of the theoretical results. 

The effective approaches emphasize the collective ori- 
gin of the plasma properties which can be best under- 
stood within a semi-classical framework. Indeed, in a 
recent work [2^ , the successes of hydrodynamics in inter- 
preting and understanding the experimental observations 
from RHIC has been reviewed. Since more exciting and 
discerning data is expected from LHC experiments soon, 
and given the above context, it is worthwhile exploring 
semi-classical techniques to understand the properties of 
QGP in heavy ion collisions. In this context, it is known 
by now [2^, [la, [13, [11] that a classical behavior emerges 
naturally when one considers hard thermal loop(HTL) 
contributions. A local formulation of HTL effective ac- 
tion has been obtained by Blaizot and lancu who have 
succeeded in rewriting the HTL effective theory as a ki- 
netic theory with a Vlasov term [1^, [s^, [3l|, Is^l . A signif- 
icant development in this direction is the realization that 
the HTL effects are, in fact, essentially classical and that 
they are much easier to handle within the frame work of 
classical transport equations ^33.,,34. ,35j. Thus, the semi- 
classical techniques appear hold the promise of providing 
tools to understand the bulk properties of QGP. 

The present paper continues the theme, and its central 
aim is to combine the kinetic equation approach which 
yields the transport properties, with the hot QCD equa- 
tions of state to make predictions which can be perhaps 
tested in heavy ion collisions. Recently, Ranjan and Rav- 
ishankar have developed a systematic approach to deter- 
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mine fully the response functions of QGP, with a special 
emp hasis on the color charge as a dynamical variable 
[iJ] . In parallel, Chandra, Kumar and Ravishankar have 
succeeded in adapting two hot QCD EOS to make pre- 
dictions for heavy ion collisions [4l| . They have shown 
that the interaction effects which modify the equations 
of state can be expressed by absorbing them into effec- 
tive fugacitics (zg,g) of otherwise free or weakly intcract- 
ing quasi quarks and gluons. Since the analysis in Ref. 
([Ijl) was illustrated only for (the academically interest- 
ing) case of ideal quarks and gluons, it is but natural to 
bring the two studies together and explore what the hot 
QCD EOS have to predict for heavy ion collisions. We 
take up this program in this paper. 

The main result of this paper is the determination of 
the modification that the heavy quark potential under- 
goes in a medium constituted by interacting QGP, as 
predicted by the two EOS which we consider. After de- 
termining the screening length as a function of tempera- 
ture, wc focus on the Cornell potentialjs^ and study the 
dissociation mechanism for cc and bb states. The results 
are rather surprising and may as well signal the inapplica- 
bility of these EOS to describe the deconfined phase. On 
the other hand, if the transition from the confined to the 
deconfined state is not a phase transition as several stud- 
ies predict [s^] , it may still be possible to attribute some 
physical significance to the predictions of these EOS. We 
undertake the project here. We show that, by using one 
of the phcnomcnological EOS is quite a good approxi- 
mation to the more rigorous lattice results, the value of 
the phenomenological coupling constant that occurs in 
the Boltzmann equation can be fixed. Ultimately, the 
physical viability or otherwise of the results need to be 
established by comparing them repeating the analysis of 
[4H with the lattice EOS. That will be taken up in a 
separate paper. 

We consider two specific hot QCD equations of state: 
The first, which we call EOSl is perturbative. with con- 
tributions up to 0((7^)[38, 39]. The second EOS has a free 
parameter 6, and is evaluated upto 0[a^ logfl/q)][4^. 
We denote it by E0S(5. d may be fine tuned to get a 
reasonably good agreement [43| with the lattice results 
[40| . which we exploit here. Both the EOS are expected 
to be vahd for T>2Tc^, and E0S6 is reliable beyond 
T ~ 4Tc. 

The paper is organized as follows: In section II, we 
introduce the two hot QCD equations of state and out- 
line the recently developed method[4l| to adapt them for 
making definite predictions for QGP at RHIC and the 
forthcoming experiments at LHC. In section III, wc ob- 
tain the expressions for the response functions of interact- 
ing QGP and in section IV, we study their temperature 
dependence in detail. In section V, we study the mod- 
ifications in heavy quark potential due to the hot QCD 
medium. We further study the temperature dependence 
of the Debye screening lengths in hot QCD. We investi- 
gate the "melting phenomena" of heavy quarkonia such 
as J/^ and bb in the medium, and extract the dissoci- 



ation temperature. In doing so we also relate the phe- 
nomenological charge that occurs in the transport equa- 
tion to lattice and experimental observables. We con- 
clude the paper in section VI. 



II. HOT QCD EQUATIONS OF STATE AND 
THEIR QUASI-PARTICLE DESCRIPTION 

There are various equations of state proposed for QGP 
at RHIC. These include non-perturbative lattice EOS 
[40| . hard thermal loop(IITL) resumed EOSjSand per- 
turbative hot QCD equations of state [H, |39|.|42|| . In the 
present paper, we seek to determine the chromo-electric 
response functions for QGP by employing two EOS: (i) 
the fully perturbative 0{g'^) hot QCD EOS proposed by 
Arnold and Zhaifll] and Zhai and Kastening [3^, and 
(ii) The EOS of 0[g^{\n{l/g) + S)] determined by Ka- 
jantie et al [4^ . by incorporating contributions from non- 
perturbative scales, gT and g^T. We employ the method 
recently formulated by Ranjan and Ravishankar [ij to 
extract the chromo-electric permittivities of the medium. 
EOSl reads 



45/34 



(1 



21^/ 
32 



15 



(1 



5jV/,as 

TT 



12 



30(1 



D TT 



(237.2 4- 15.97iV/ - 0.413iV 



135 



/ 



(1 



x(l + ^))_(l 
^ 6 '^'^ 8 ^ 



12 



33 



)ln[ 





ln(^ 

TT 




TT 


27r ^ 



799.2- 21.99iV/ - 1.9267V/ 



495 , Nf-, 
+OiasfHas)). 
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As mentioned earlier, S is an empirical parameter, in- 
troduced to incorporate phcnomcnologically the undeter- 
mined contributions at 0{g^). It also acts as a fitting 
parameter to get the best agreement with the lattice re- 
sults. . 
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FIG. 1: (Color online) Relative equation of state(uiri ideal 
EOS) for pure gauge theory plasma as a function of T/Tc for 
various values of S. 
0.9 




FIG. 2: (Color online) Relative equation of statew;rt ideal EOS 
for full QCD plasma with Nf = 2, 3 as a function of T/Tc for 
various values of 5. 



A. The underlying distribution functions 

The construction of the distribution functions that un- 
derUe the EOS, in terms of effective quarks and gluons 
which act as quasi-excitations, has been discussed by 
Chandra et. al., [4l[ in the specific context of EOSl and 
E0Si5. To review the method briefly, all the terms that 
represent interactions are collected together by recasting 
them as effective fugacities {zq^g = exp(^q g)) for the oth- 
erwise free quarks and gluons. Of course, the pure gague 
theory case is simply obtained by putting the number 
of flavors, Np = in the EOS. Thus, /ig represents the 
self interactions of the gluons, while ^/ encapsulates the 
quark-quark and the quark gluon interaction terms. Im- 
portantly, the two EOS of interest to us are valid when 
T > 2Tc, and in this range, the quantities = P^q.g 
are perturbative parameters. Thus, it is possible to solve 
for jif\g self consistently through a systematic iterative 
procedure. In this procedure, all the temperature effects 
are contained in the effective fugacities z = z{as{T/Tc)), 
where we display the dependence on the temperature and 
coupling constant explicitly. It has been shown in Ref. 
[m (where the details can be found) that one can trade 
off the dependence of the effective fugacities on the renor- 
malization scale ( fJ-Ms) by their dependence on the crit- 



ical temperature Tc- For that purpose, one utilizes the 
one loop expression of as(T) at finite temperature given 
by [13 
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employing which the dependence on ^7^75 is eliminated, 
in favor of Tc- Consequently, the effective chemical po- 
tentials get to depend only on T/Tc- Note that effective 
fugacities have merely been introduced to capture the in- 
teraction effects present in hot QCD equations of state. 

Once the distribution functions are in hand, the study 
of transport properties is a straight forward exercise if 
we employ the analysis put forth by Ranjan et. al. [l^ . 

In Figs. 1 and 2, we display the behavior of EOSJ for 
various values of the parameter S. The figures show the 
pure gauge theory contributions to the EOS and full QCD 
separately. We remark parenthetically that the studies 
in the earlier work [4l| were confined to EOSl and the 
special case (5 = in E0Si5. For the details on EOSl and 
E0S(5 for (5 = 0, we refer the reader to Ref. [4l[ (see Fig.l- 
7 of Ref.jlll]). First of all, we see that as 6 increases in 
magnitude, the EOS, for both pure gauge theory and full 
QCD, become softer, with P/Pi taking smaller values, we 
denote the the ratio P/Pj by Kajantieji^l obtains 
the best fit with the lattice results of Boyd et. al.jil] by 
choosing a value S ~ 0.7. We find that to get agreement 
with the more recent results of Karsch [40], S w 1.0 is 
preferred, when we consider T > 2Tc. In short, we find 
that the range of values 0.8 < S < 1.2 gives a reasonably 
good qualitative agreement with the lattice results for 
the screening lengths. 

Here, we wish to mention that there is an uncertainty 
in fixing the free parameter 6. This follows from the 
freedom in choosing the QCD renormalization scale at 
high temperature. This has been investigated in detail 
by Blaizot, lancu and Rebhan 45[. The value of 5 in the 
present paper has been obtained by employing the one 
loop expression for the running coupling constant and the 
QCD renormalization scale determined in Ref. [l^. We 
intend to study the qu asi-p article content of HTL and 
HDL equations of state [U Ell and lattice equation of 
state in future. 

The behavior of the corresponding fugacities, as a func- 
tion of temperature, is shown in Fig. 3. It may be seen 
that < Zg^q < 1.0 which ensures the convergence of the 
method to determine the effective fugacities from the hot 
QCD EOS. We now proceed to determine the response 
of the plasma in the next section. 
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FIG. 3: (Color online) Efltective parton fugacities (zg,,) quarks 
determined from E0S(5 as a function of temperature. Note 
that the behavior is shown for S = 1.0. 



III. RESPONSE FUNCTIONS FOR 
INTERACTING QGP 

Recently Ranjan and Ravishankar [ij] have deter- 
mined the form of chromo-electric response functions for 
collision less quark-gluon plasma within the framework 
of semi-classical transport theory. They have set up the 
transport equation in the extended phase space includ- 
ing the SU(3) group space corresponding to dynamical 
color degree of freedom. They have taken the distri- 
bution function in a coherent state basis defined over 
the extended single particle phase space TZ^ (8) Cg , where 
Cg = GfH is the phase space corresponding to the color 
degree of freedom, obtained as a coset space by factoring 
the group space by the stabilizer group H of any refer- 
ence state in the Hilbert space. Having been employed to 
study the ideal case, the formalism has not been applied 
to examine the behavior of the plasma with a realistic 
EOS. We employ the results of the previous section and 
rectify this drawback, by incorporating the interaction 
effects as represented by EOSl and E0S(5. 

A brief comment on the response functions. In contrast 
to electrodynamic plasma, the chromo-electric response 
has a richer structure. Apart from the standard permit- 
tivity which we shall call Abelian and denote by e^, there 
are additional response functions, their number depend- 
ing on the color carried by the partons. Thus, quarks 
have an additional response function which affects the 
non- Abelian coupling. The corresponding permittivity 
will be called non- Abelian, and denoted by e^r. The two 
functions exhaust the response in the quark sector. The 
gluonic sector, arising from the adjoint representation of 
the gauge group admits yet another kind of response, 
corresponding to tensor excitations. These excitations 
are not allowed in the quark sector (which emerges from 
the fundamental representation of the gauge group). We 
consider each of these response functions for the inter- 
acting QGP. The response functions are obtained in the 
temporal gauge. 

Consider first the familiar Abelian component of the 
response e^. For an isotropic plasma(in the absence of 



chromo-magnetic fields), its expression is given by 

~eAiiu,k) = l + Q^Ioiiu,k) (4) 

where = Q"^Q'^ is the color charge magnitude squared, 
and Iq is determined by the equilibrium distribution func- 
tion thus: 
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The non-Abclian response function, which has been 
evaluated in the long wavelength limit, is given by 



e^iuj^uj') = {1 + 
where Ii is defined as 

h{oj,k)^^Tr 
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We recall that the new constitutive Yang-Mills equations, 
in the presence of the medium, arc given by 
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As pointed out in the Abelian and non- Abelian 
responses are not independent of each other. Gauge in- 
variance relates them, by virtue of which we can obtain 
both from a common generating function as follows: 



1 d f k ■ p 
^o = -j^-Q^ / '"'■'^ —)hdpjeqd^p 

= 'l^'iA] I ^"^"^ ~ ^)dpJe,d'p^ ■ (8) 

We further recall that these expansions are determined 
when the system is displaced slightly from its equilib- 
rium, in the coUisionlcss limit. 



A. Ideal response 

It is convenient to first write the expressions for the re- 
sponses of ideal distributions for quarks and gluons. The 
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responses due to EOSl and E0S5 get a simple modifica- 
tion over their ideal forms since we have mapped success- 
fully the interaction effects into quasi free partons with 
effective fugacities. Thus, in the ideal case we have, for 
the quarks. 



2}] 
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and the non-Abelian response function is given by 



= {1 - 
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The imaginary part of Abelian(e^) and non-Abelian 
component (eAc) of the chromo-electric permittivity can 
be easily evaluated by the standard Landau ie prescrip- 
tion. These are needed to obtain landau damping which 
we do not study here. 

The contribution to the permittivity from the gluons is 
closely related, and not independent of the contribution 
of the quarks written above. Indeed, if we define the 
susceptibilities 



and 



1 



for the quarks and the gluons. It can be shown that [lj| 
the gluonic permittivity can be simply read off from the 
quark permittivity (and vice versa) as 
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where Np is the number of flavors. In short, for the 
total susceptibility, we have the simple relation Xq'^ — 

Np A.N 
2 ' 
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= 

via the integral below 
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The function g'2{zg) = ^g2{zg) where g2{zg) is defined 
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gu{zg) has the series expansion 
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Note that 52(1) = 1 gives the ideal limit. 

Similarly, the corresponding expressions for in the 
quark sector are obtained as 
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and the non-Abelian response for effective quarks 
reads: 



^^^Q^T^Nff^{zf) 1 
ejv = {1 71- 
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(15) 



The function /2(^/) = ^f2{zf) where f2{zf) is defined 
via the integral below. 
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B. Interaction effects 



and /^(l) = 1. 



We now consider the modification that the above ex- 
pressions undergo permittivities arising because of the 
new EOS. Recall that the corresponding equilibrium dis- 
tribution functions differ from each other only in their 
form for the chemical potentials Pq.g. The responses thus 
depend on the interactions implicitly through an explicit 
dependence on Zq,g. 

Considering the gluonic case, i. c, pure gauge theory 
first, we get the expressions for the two permittivities as 



w - 
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and the non-Abelian response function is 



IV. EFFECTIVE CHARGES AND RELATIVE 
SUSCEPTIBILITIES 

Eq. (|12j|l3|) admit a simple physical interpretation, 
when compared with their counterparts Eg . ([^ - fTU|) . In- 
deed, the sole effect of the interactions on the transport 
properties is to merely renormalize the the quark and the 
gluon charges Qg^q as shown below: 



'92{zq); 



The renormalization factors g^^Zg)., f2i^f) further possess 
the significance of chromo-electric susceptibilities, rela- 
tive to the ideal values. To see that, we note that the 
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Abelian and the non-Abelian strengths for gluons as well 
as quarks suffer the same renormalization reflecting the 
underlying gauge invariance. Furthermore, the expres- 
sions for the relative susceptibilities are given by, 



7^ 



_ X{z) ^ A{z) _ U{z) _ ( f!^{zf) for quarks 



X(l) ^(1) 1.92(^9) for gluons 



(16) 



and 



7^o 



(17) 



Note that the relative susceptibilities are entirely func- 
tions of the single variable T jT^-, and are independent of 
{ijj,k). The dependence of the susceptibilities on {u!,k) 
has already been studied in detail in Rcf. . We merely 
concentrate on the temperature dependence below. 

Before we go on to discuss the susceptibilities and other 
bulk properties, we point out an essential care to be taken 
in using the above susceptibilities for determining the 
response of the plasma. For pure gauge theory, only the 
giuonic part contributes, while for the full QCD, we have 
to necessarily take the contribution from both the quark 
and the giuonic sector. We discuss both the cases below. 
The response functions for the full QCD is obtained by 
averaging up the above calculated response functions for 
quark as well as gluon plasma. The relative susceptibility 
for full QCD plasma is given by 



7^' 



- ^(1) 

^ Nff^{zf) + 2g'^{z,) 



(18) 



where z is the effective fugacity of partons in full QCD 
plasma. 



A. Behavior of the susceptibiUties 

We now proceed to study the behavior of the relative 
susceptibilities displayed in Egs. pB]) . ([T7]) and (fT5|) as 
functions of temperature. As observed, relative suscepti- 
bilities for both quarks and gluons scale with T/Tc. We 
have plotted the relative susceptibilities TZ, TZqg and TZ' 
as functions of T/T^ (See Figs.4-7), for both EOSl and 
E0S(5. Please note that we have chosen S = 1.0 in E0S(5. 

Fig. 4 shows the relative susceptibility of a purely giu- 
onic plasma as a function of temperature for EOSl and 
E0S(5. 

We see From Fig. 4 that the susceptibility of a purely 
giuonic plasma is weaker in the presence of interactions, 
approaching its ideal value asymptotically with increas- 
ing temperatures. Equivalently, there is a decrease in the 
value of the phenomenological coupling , relative to its 
ideal value. 

The behavior of quark gluon plasma is not qualitatively 
different from that of a purely giuonic plasma, as may be 




FIG. 4: (Color online) Relative susceptibility, g'2{zg) (see Eq. 
(|16() . for pure gauge theory plasma as a function of T/Tc for 
EOSl and BOSS [5 = 1). 




FIG. 5; (Color online) Relative susceptibility, defined in Eq. 
(HSl), for the full QCD plasma as a function of T/Tc, for EOSl 
and EOS(5 {S — 1). We have studied the cases Nf = 2, 3. 



seen from Fig. 5. In other words, the quark contribu- 
tion is of the same order as the purely giuonic contribu- 
tion. However, the relative contribution from the quarks 
and the gluons does depend on the EOS considered. In- 
deed, with EOSl (where interactions up to 0{g^) are 
included), Fig. 6 shows that the quark contribution dom- 
inates slightly over the giuonic contribution for Np ~ 2. 
The dominance is more pronounced for the more realistic 
case Np = 3. In contrast, we sec from Fig. 7, that E0S(5 
(with (5 = 1) predicts that the giuonic contribution is 
marginally larger for Np = 2 and becomes sub dominant 
when Np = 3. This distinction between the two EOS is 
of no practical consequence since, given Tc ~ 170MeV, 
one has to necessarily work with Np = 3 at T = 2Tc. . 




FIG. 6: (Color online) Ratio of the quark to gluonic contribu- 
tions to the susceptibility (see Ea.(|17p as a function of T/Tc, 
as predicted by EOSl. 
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FIG. 7; (Color online) Ratio of the quark to gluonic contribu- 
tions to the susceptibility (see Ea.(|17p as a function of T/Tc, 
for EOS<5, 5 = 1.. 



V. THE HEAVY QUARK POTENTIAL 

Now we shall apply the results of the previous sections 
to discuss the heavy quark potential in the presence of 
interacting medium. We consider the Cornell potential 
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At 



where a and A are phenomenological constants. The 
first term shows the Coulombic behavior and dominates 
at small distance while the second term causes linear con- 
finement, dominating at large distances. 

It had been expected earher that the long range part 
of the Cornell potential does not survive in the quark 
gluon phase. This expectation assumes a phase transi- 
tion from the hadronic to deconfined phase. More recent 
studiesjlj indicate that in all likelihood, deconfinement 
is not a phase transition, but a crossover. If such to be 
the case, there is no reason to expect the linear part of 
the potential to disappear completely. With this in mind, 
we study the modifications of both the Coulomb and lin- 
ear terms, and examine how reasonable the EOS under 
consideration are. 



Since the potential has no explicit color dependence, it 
is sufficient to employ the Abelian components of the per- 
mittivities. At a; = 0, the quark and gluon permittivities 
have the form 



e,(fc,T) 
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Therefore the full permittivity reads 



= 1 

= 1 
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in terms of the Dcbyc mass ?ti|) = SttQ^T'^ [iV//2(zq) -|- 

92{Zg)\. _ 

The qq potential undergoes a modification due to the 
medium via e{k,T), as given by 4'{k) — > (j){k)/e{k,T) = 
(j)s{k,T). We note that in determining the Fourier trans- 
form of Cornell potential, we regulate the linear term 
exactly the same way we regulate the Coulomb term, by 
multiplying with an exponential damping factor. The 
damping is switched off after the Fourier transform is 
evaluated. The Fourier transform is thus obtained as 
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The modified potential thus acquires the form 
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(21) 



(22) 



We note that for a gluonic plasma, m|j = 
167:Q^T^g2{zg). 

On comparing Eq[22]with Eq[2T]we infer the rcnormal- 
ization of the couplings 



A, 



A 



fc2 



^eff 



fe2 



A. Screening of the heavy quark potential 

Of interest to us is the form of the potential in the real 
space, as a function of spatial separation. The inverse 
Fourier transform yields it to be 



2A 
ml 
2A 



a) 



exp {^mjjr) 



2A 

rriD 



(23) 
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It follows from the above equation that the medium 
transforms the linear potential to the long range Coulomb 
form, just as it modifies the bare Coulomb term to the 
short ranged Yukawa.. The modified potential is not 
short ranged; it is not confining cither. To appreciate 
this, note that at large T, the above expression reduces 
to 



2A 

— 2 am£i 



(24) 



Thus, contrary to the MaxwcUian plasmas which support 
only short range interactions, the two EOS predict that 
the heavy quark potential continues to be long ranged, al- 
though absolute confinement, which was a quintessential 
feature of the unscreened potential, is lost. It might as 
well be that the above results signify that the EOS fail to 
describe the hadronic matter in its deconfined state. On 
the other hand, since the transition from the hadronic to 
QGP phase could be a cross over, and not a phase tran- 
sition [33, it could be possible that the above result is 
not entirely devoid of physical significance. If we adopt 
the latter view, if only for the purposes of analysis, a 
discussion of screening cannot, therefore rely entirely on 
the interpretation of inverse Debye mass as a screening 
length in its usual sense. We address the issue below . 

Let us consider the high temperature limit of the po- 
tential, given by Eq. ([M)) . Ignoring the additive contribu- 
tion, the energy of the qq in the ground state is simply 
given by 



— J— 



where niq is the mass of heavy quark. 




FIG 
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. 8: (Color online) Debye screening length for gluonic 
quark-gluon plasmas as a function of T/Tc for EOSl. 
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FIG. 9: (Color online) Debye screening length for pure gauge 
theory plasma in EOS(5 as a function of T/Tc for various values 
of 5. 



by the matching conditions. In the case of pure gauge 
theory. 



^ = (5^)^5l(%). 



(25) 



3847r2Q4y5 ' T, 
And for fuU QCD: 

^£$fl = li^n^fM^.)+9.i^.)] (26) 



B. Estimation of Q and a determination of the 
screening length 

The above equation is still not amenable to comparison 
with experiments since it has the undetermined parame- 
ter Q. To estimate Q, we need an additional input which 
we obtain by comparing the screening length obtained as 
a solution of Eqs. (pS)) and with the lattice results, 
reported by Kazmarek and Zantow ji^. Note that the 
screening lengths for gluonic and quark gluonic plasmas 
are respectively given by 



\s — 



1 



1 



QT 



An = 



16'Kg2{zg 



1 



(27) 



QT 



i^87riNfh{z,)+g2{z,)) 



(28) 



The binding energy is, of course, temperature depen- 
dent and approaches zero as T ^ 00. At any finite tem- 
perature though, the quarks possess a thermal energy 
Erp ~ T (by equipartition theorem), leading to an ion- 
ization of the quarkonium when Ex matches the binding 
energy. The dissociation temperature is determined 



Recalling that our results are valid for T > 2Tc, we 
match the pure gauge theory result with the lattice values 
A - 0.15 fm, and Tc = 0.27 GeV. We obtain 



0.19 



1 



1 



o.27 4g(r/re) y(7rg2(^s))' 



(29) 
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This leads to the estimate Q ~ 0.15. The temperature 
dependence of the screening lengths can be thereafter 
determined for the two equations of state. We emphasize 
that the choice 6 = 0.9 gives the best agreement between 
E0Si5 and the lattice EOS for gluonic plasma [ioj . 



C. Dissociation temperatures for quarkonia 

Since there are no free parameters left, it is a straight 
forward task to determine the dissociation temperatures 
for the heavy quark bound states. We are principally in- 
terested in J/'^ and bb states, for which we have gathered 
the results in Table. 1, after obtaining graphical solutions 
for Eq. ([^ and Eq.f^H). We have employed the values 
iTLc = 1.5GeV, rub = A.^GeV and A = O.lSGe^^ for the 
quark masses and the strength of the Cornell potential. 
It is noteworthy that the dissociation temperatures are 
all roughly in the range Td w (2 — 3)Tc, which is higher 
than the temperatures achieved so far. Since the temper- 
atures expected at LHC is in the range T ^ 2Tc — 3Tc, 
one may expect to test these predictions there. 

We now turn our attention to compare hot QCD es- 
timates for dissociation temperatures with other theo- 
retical works. In a recent paper, SatzjS^ has studied 
the dissociation of quarkonia states by studying their in- 
medium behavior. These estimates were based on the 
Schrodinger equation for Cornell potential. In a more 
recent work, Alberico et aZ[4^ reported the dissociation 
temperatures for charmonium and bottomonium states 
for Nf ^ and Nf = 2 QCD. In this work, they have 
solved the Schrodinger equation for the charmonium and 
bottomonium states at finite temperature in the presence 
of temperature dependent potential- computed from the 
lattice QCD. We have quoted these results in Table 2. 
The estimates for Nf = and Nf = 2 cases for both 
EOSl and E0S<5(Table 1) arc closer to Ref.IH. On the 
other hand the estimates for J/\l/ dissociation temper- 
atures for both EOSl and E0S2 are larger than that 
of Ref. (49| while bottomonium dissociation temperature 
estimates are slightly smaller. We do not have lattice 
estimates at present to compare the dissociation temper- 
atures for Nf = 3 QCD. However the hot QCD esti- 
mates are consistent with the lattice predictions [50| on 
the survival of heavy quarkonia states near 2Tc and pre- 
dictions of dynamical Nf — 2 QCD by Aarts et a^|5ll|. 
Along these results, we wish to mention the very recent 
estimates on dissociation temperature reported by Mocsy 
and Petreeczkyfs^l- Their estimates for J/^* dissociation 
temperature is 1.2rc and for T is 2Tc. The estimates for 
both EOSl and E0S(5 are larger as compare to these re- 
sults. 



1. Comparison of Debye screening length with lattice results 

Finally, with a view to benchmark our estimates of the 
screening lengths, by comparing them with the recent lat- 



TABLE I: The dissociation temperature(rD) for various 
quarkonia states (in unit of Tc). 



Hot EOS 


Quarkonium 


Pure QCD 


iVf = 2 






T /\lr 

J/ y 










T 


9 ^ 






EOS5 


J/* 


1.86 


2.38 


2.24 


(5 = 0.8 


T 


2.12 


2.76 


2.58 


EOS<5 


J/* 


1.95 


2.45 


2.32 


5 = 1.0 


T 


2.2 


2.83 


2.66 


EOS5 


J/* 


2.03 


2.52 


2.40 


5 = 1.2 


T 


2.28 


2.9 


2.74 



TABLE II; The dissociation temperaturefTp ) for various 
quarkonia states (in unit of Tc) from Ref.^^ and Ref. [49l| . 
The first and third rows are the estimated values for disso- 
ciation temperature from Ref. [i^ and second and fourth are 
from Ref. [49| 



Quarkonium 


^/=0 


Nf = 2 


J/* 


2.1 


>2 


J/* 


1.40 


1.45 


T 


>2 




T 


2.96 


3.9 



tice results reported by Kazmarek and Zantow |47| . we 
plot 2Xd as a function of T/Tc- The results are shown for 
EOSl as weU as E0S,5, for three values 5 = 0.8, 1.0, 1.2. 
The results for pure gauge theory (gluonic plasma) are 
shown in Fig. 9, and Fig. 10 shows the results for full QCD. 
We find that on comparison with Fig. 2 of Ref. [i^l, these 
values, (5^1 are the most favored which justifies, a pos- 
teriori our choice for the parameter. 




2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 



T_ 

FIG. 10; (Color online) Debye screening length for full QCD 
plasma as predicted by E0S5 as a function of T /Tc for various 
of 5. 
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VI. CONCLUSIONS AND OUTLOOK 

In conclusion, we have successfully extracted the quasi- 
free particle content of two hot QCD equations of states 
and used them to determine the chromo-electric permit- 
tivities within the standard Boltzmann-Vlasov kinetic 
approach. The Abelian and the non-Abelian components 
of the permittivity are obtained, for pure gauge theory 
and the full QCD. We have shown that the effect of the 
interactions is to merely renormalize the magnitude of 
the effective color charge, Q. We have used the permit- 
tivities to study critically the modifications in a realistic 
heavy quark potential. The dissociation temperatures 
are carefully estimated, by fixing the magnitude of Q by 
an explicit matching with a lattice result. The values 
obtained are quite close to the exact lattice results. The 
viability of the two EOS, especially E0Si5 is thus phe- 
nomenologically well supported. Our analysis suggests 



strongly, and in agreement with the lattice results, that 
J/'f suppression can be seen in QGP only for T > 3Tc. 

A true test of the above predictions would be possible if 
we succeed in extracting a quasi particle description from 
the lattice EOS. Studies are under way in this direction. 
It should also be of interest to extend the analysis to 
other signatures like strangeness enhancement, and also 
for QGP with a finite baryonic chemical potential (ssl . 
[S^ l . These will be taken up in a later work. 
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